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1 Introduction 

Let J\f(x) denote the number of positive integers n < x for which n 2 + 1 is 
square-free. It was shown in 1931 by Estermann [I] that 



for x > 2, where 



N{x) = c x + 0(x 2/3 log x) 



C0 = 2 

p=l mod 4 



Estermann's argument is very simple, but despite the passage of 80 years 
the exponent 2/3 appearing above has never been improved. The aim of the 
present paper is to establish the following result. 

Theorem We have 

N(x) = c x + £ (x 7/12+£ ) 

for any fixed e > . 

It is easy to construct intervals (x, x + c log x] with a small positive con- 
stant c, such that n 2 + 1 has a non-trivial square factor for every n in the 
interval. This shows that the error term in our theorem is Q(\ogx). However 
we know of no better result of this type, and it is unclear what one should 
conjecture. With the much simpler problem of the number of square- free 
integers n < x one has an easy error term 0{x 1 ^ 2 )^ but any reduction in the 
exponent 1/2 would appear to require a quasi Riemann Hypothesis. Thus it 
seems unlikely that we could reduce the exponent in our theorem below 1/2 
without a radically new idea. 
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The key point in our treatment will be to give good upper bounds for the 
frequency of solutions to the Diophantine equation d 2 e = n 2 + 1. Analysing 
this over Q(i) we are led to study the condition 2x1^22/1 + (%\ — x\)y2 = 1, 
which we may interpret as saying that the point (s,t) = (^1/^2,2/1/2/2) lies 
close to the curve t = (1 — s 2 )/{2s). In order to study this we we will use 
a variant of the "Determinant Method" , developed from the author's papers 



The author was introduced to this problem by Dr Tim Browning. His 
contributions to the resulting discussions, and his careful proof-reading of 
the present paper, are gratefully acknowledged. 

2 Preliminaries 

For the proof it will clearly suffice to show that 



The argument begins by observing that for x > 1 and 1 < D < x we have 



0, ®. 



M{2x) - N{x) = c x + £ {x 7/12+£ ) 



M{2x) - M{x) 



x<n<2x 



E E ^ 



x<n<2x d 2 \n 2 +l 



E < n < 2x : d 2 \ n 2 + 1} 




d<D 



+ Q{ E #{x<n<2x:d 2 \n 2 + 1}). (1) 



D<d<4x 



For d < D we write 



p{d) 



p = #{m mod d 2 : d 2 \ m 2 + 1} 



2 



and we take mi, ... , m p to be a corresponding set of admissible values for m. 
Then 

p 

#{x < n < 2x : d 2 \ n 2 + 1} = ^ #{x < n < 2a; : n = mj mod d 2 } 

= x^- + 0(p(d)). (2) 
Thus terms with d < D contribute to (1) a total 

xJ2v(d)p(d)d- 2 + 0(J2p(d)). 

d<D d<D 

The function p(d) is multiplicative, with p{p) = 2 for p = 1 mod 4 and 
pG°) — otherwise. Thus p(cf) is bounded by the familiar r(rf) function 
which counts representations as sums of two squares. We therefore see that 

E<d<2E 

for any integer E, whence 

E \p{d)p{d)d- 2 \ < IT 1 and ^ < D - 

The contribution to (1) corresponding to values d < D is therefore 

x ±tm^ + ( x D-') + 0(D). 
d=l 

Since 

d=l p 

we see that this produces the main term in our theorem. We will minimize 
the other error terms by choosing D = x x l 2 . 

To handle the larger values of d we consider dyadic ranges E/2 < d < E, 
and write 

M(E, F) = #{(e, /, n) G N 3 : E/2 < e < E, F/2 < f < F, e 2 f = n 2 + 1}. 
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Then the range d > D contributes to (1) a total 

< V max M(E, F), 

where the summation for E runs over powers of 2. Thus our problem reduces 
to one of estimating Ai(E, F) efficiently. Heuristically one might expect that 
e 2 f — 1 is a square "with probability" of order (e 2 /) -1 / 2 . This leads one to 
conjecture that the true order of magnitude for M(E,F) might be about 
F 1 / 2 . For his proof, Estermann showed that 

#{(e,n) G N 2 : E/2 < e < E, e 2 f = n 2 + 1} < logx, (3) 

whence Ai(E, F) <C F logx. One easily sees how this leads to the error term 
O (x 2 / 3 logx). We will need sharper bounds, but we note that Estermann's 
estimate shows that the range F < x 1//2 yields a satisfactory contribution. 
Since we have taken D = x x l 2 we may therefore assume in what follows that 
x 1/2 < E < x 3/4 and x l/2 < F < x. 



3 The Determinant Method 

We begin our analysis of Ai(E, F) by using the unique factorization property 
for This shows that if e 2 f = n 2 + 1 then there are integers Xi,x 2 ,yi, yi 
for which 

e = xl + xj, f = yf + y 2 and (x 1 + ix 2 ) 2 {yi + iy 2 ) = n + i. 

It follows on taking the imaginary part that 2x\X 2 yi + (x\ — x\)y 2 = 1. If 
> \x 2 \ we will swap x\ and x 2 , and change the sign of y 2 . Hence we 
may suppose, without loss of generality, that < \x 2 \, and hence that 
\xi\ < E 1/2 and E l/2 < \x 2 \ < E 1 ' 2 . We observe that 

max{|2xiX2| , \x\ — x 2 ,]} ^> x\ + x\ ^> E. 

Thus if we write qi(xi, x 2 ) = 1x\X 2 or x\ — x\ as appropriate, and take q 2 to 
be the alternative quadratic form, we may assume that E ^ \qi(xi, x 2 )\ <C E 
and 52(^1,^2) ^ E. Then, labelling yi,y 2 either as -21,-22 or as z 2,Zi w e will 
have 

qi(xi,x 2 )z 1 + q 2 (x 1 ,x 2 )z 2 = l, (4) 

whence 

E\ z i\ < \qi( x i, x 2) z\ I < 1 + ^2(^1, ^2)^2! < 1 + -E 1 1 ^2 1 • 
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Since we cannot have z\ = z 2 = we deduce that \z\\ <C \z 2 \. Then, since 
F < z\ + z\ < F we see that \z x \ < F 1 ' 2 and F 1 ' 2 < \z 2 \ < F 1 ' 2 . 
We now deduce from (jl]) that if s = x\/x 2 and t = z\j z 2 then 

t = + 0(E-'F^ 2 ) = -^M + 0(a .-i). 

Thus if we write = — q 2 (s, l)/gi(s, 1) then the point (s,t) lies close to 
the curve 4>(S) = T. Our task is therefore to estimate the number of rational 
points (s,t) with s,{< 1 lying within 0(x~ l ) of the curve <fi(S) = T, and 
for which the "heights" of s and t are at most E 1 ^ 2 and F 1 ^ 2 respectively. 

The situation here is similar to that in the author's paper [6]. We shall 
use a real- variable version of the "determinant method", but there is an 
important difference, in that the variety given by the equation 

q 1 (x 1 ,x 2 )z 1 + q 2 {x ll x 2 )z 2 = 

lies naturally in P 1 x P 1 , rather than in P 2 . Indeed this makes our situation 
correspond exactly to that considered by Huxley [7] , [S] • Unfortunately Hux- 
ley's bounds, which were obtained for general plane curves, are not strong 
enough for our application. In particular, he focuses on the case in which </> 
is not a rational function. 

Following the method from the author's work [6, §2] we choose an inte- 
ger parameter M £ [x 1//2 ,x] and split the available range for s into O(M) 
subintervals / = (sq, Sq + M -1 ]. We then investigate the number of solutions 
in which s belongs to a particular interval J. If we write s = sq + u we find 
from Taylor's Theorem that <j)(s) = 0(so) + u<p'(so) + 0(x _1 ). Hence if we 
set v = t — 0(s o ) — u(f>'(s ) we will have s = s + u, t = 0(s o ) + u(f>'(s ) + v 
with v x" 1 . We now label all the solutions corresponding to the interval 
I as (sx, ti), . . . , (sj, tj), say. We proceed to choose positive integers K, L 
and to label the monomials s k t l for k < K, I < L as mi(s,t), . . . ,rriH{s,t), 
where H = (K + 1)(L + 1). The determinant method uses the J x H matrix 
M, whose jh entry is rrih(sj,tj). The aim is to show that the rank of Ai is 
strictly less than H. If this can be achieved, one may deduce that there is a 
non-zero vector c with 

Mc = 0. (5) 

This vector c may be constructed out of appropriate subdeterminants of M.. 
Thus its entries will be rational numbers with numerators and denominators 
of size <^k,l x H ( K+L ^ , since Sj and tj have numerators and denominators of 
size x 1 / 2 x. We now observe that the matrix equation §5§ means that 
there is a polynomial C(s, t), with coefficients given by the vector c, such that 
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C(sj,tj) = for all pairs (sj,tj). Multiplying out the common denominator 
of the coefficients we may assume that C has integer coefficients, of size 

This is one of the key stages in the proof. We deduce that all points (s, t) 
for which t is close to <f>(s), and for which s lies in an appropriate short range 
I, actually lie on the curve C(s,t) = 0. 

We now show that M. does indeed have rank less than H, if the parameter 
M is suitable chosen. For this we select any H x H sub determinant, A say, 
from Ai, and show that A = 0. Without loss of generality we may suppose 
that A comes from the first H rows of Ai. Since the j-th row contains 
rationals with a common denominator of x^jZ^j it is clear that 



J<H J 



\J1 

Thus to show that A = it will suffice to prove that 

A E- KH ' 2 F- LH I 2 (6) 

with a suitably small implied constant. 

When we substitute s = so + u and t = 0(so) + u(p'(so) +v the monomials 
rrij(s, t) produce polynomials in u, v. Thus A is a generalized van der Monde 
determinant. If Uj and Vj correspond to Sj and tj then we have \uj\ < M~ l 
and \vj\ < V^ 1 for some V of exact order x. An estimate for the size of 
A is now provided by Lemma 3 from the author's work [5]. If we order all 
possible monomials M~ k V~ l in decreasing size as 1 = Mo, Mi, . . . then the 
lemma shows that 

H 

a <€. H n M h- 

h=l 

If M H = W~ l then M~ k V~ l > M H if and only if 

HogM + /logy < log W. (7) 
The number of such pairs k, I is 

vo g wr +0 /iog^' 1 + 



2 (log M) (log V) \\ogx 
and since this must equal H we deduce that 



log W = H 1/2 v / 2(logM)(\og V) + 0(log 



x 
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Moreover 

H 



\og\\M h = -^(fclogM + l\ogV) 



h=l kl 



(\ogWf + Q ((\ogWf 



3(logM)(logF) V lo & x 
the sum over k, I being subject to d7j). It follows from (jSJ) that 

log TT M h = -H 3 ' 2 ^ v/(logM)(log^) + O (# log x) , 

and hence that 

log | A| < H (I) ~ # 3/2 ^V(logM)(log V) + 0(H\ogx). 
This will be sufficient for fl6]) providing that 

y log E + | log F < (KL) x ' 2 ^ V(logM)(logy) + O k , £ ( 1 ) + O (log x) . 

In order to use this optimally we will take K = [L(log F) / (log E)]. Since our 
size constraints on E and F imply that L <C K <C L it then suffices that 

LlogF < L^v/(logM)(logy)^p5 + O l (1) + 0(\ogx). 
3 Vlog £ 

Hence if 5 > is a small positive constant, and 

^ V(logM)(logV0^P5 > ( 1 + 5) log F 
3 vlog 

it will be enough to have L = L(5) sufficiently large, and x 3><5 1. The 
condition may be rewritten in the form 

logM >i(i + ^»i£) 

8 log V 

and since K>iwe may summarize our conclusions as follows. 
Lemma 1 Let rj > be given, and suppose M 6 [x 1//2 ,x] satisfies 

(log E) (log F) 



\ogM> 1(1 + 7]) 



logx 



Then for any interval I = [s , s + M *] there is a corresponding non-zero 
integer polynomial Cj(s, t) satisfying 

C I {x 1 /x 2 ,z 1 /z 2 ) = (9) 

for any solution of Op with x\jx 2 £ /. Moreover Cj has total degree V (1), 
and coefficients of size v (x K ) for some constant k = k(t)). 

4 Counting Solutions of Equations 

While the previous section involved the application of a general method, the 
next stage in the proof requires an ad hoc argument, to count points which 
simultaneously satisfy both fll]) and (EJ). We begin by showing that it suffices 
to assume that Cj is absolutely irreducible. Let (s,t) be a rational point 
satisfying F(s, t) = for some monic factor F of Cj which is not defined 
over Q. Then F a (s,t) = for every conjugate F a . The number of possible 
points s,t is then O v (l) by Bezout's Theorem. Since X\,x 2 are coprime, and 
similarly for Zx,z% we obtain O v (l) solutions this way. Thus we need only 
consider absolutely irreducible factors F of Cj which are defined over Z. The 
height of any such factor is again bounded by a power of x, by Gelfond's 
Lemma (see Bombieri and Gubler [TJ Lemma 1.6.11] for example). Moreover 
the number of different factors to consider is O v (l). Thus it suffices to con- 
sider the case in which F{x\/xi,z\/ Z2) = for some absolutely irreducible 
polynomial F satisfying the same conditions as Cj. 

Our next move is to clear the denominators x 2 and z 2 so as to replace the 
equation F(s,t) = by a bi-homogeneous one 

F(x 1 ,x 2 ;z 1 ,z 2 ) = 0, (10) 

say. For a given interval / we will have 

s < s = — < s + M~ l . 

x 2 

It therefore follows that \x\ — sqx 2 \ < E X I 2 M~ X , since \x 2 \ < E 1 ^ 2 . If we let 

A = {(E- 1/2 M( Xl - s x 2 ) , E~ 1/2 x 2 ) : (x u x 2 ) £ Z 2 } 

then A is a lattice of determinant E~ X M, and we are interested in points 
(«i, a 2 ) £ A falling in the square 

S = {(ai,a 2 ) : max(|ai|, \a 2 \) < 1}. 
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Let g*- 1 ** be the shortest non-zero vector in the lattice and the shortest 
vector not parallel to g . These vectors will form a basis for A. Moreover we 

have Aig (1) + A 2 g (2) E S only when |Ai| < |g (1) r x and |A 2 | < |g (2) | _1 . These 
constraints may be written in the form |Aj| < Lj, for appropriate bounds 
L X ,L 2 . Since |g (2) | > |g (1) | and |g (1) |.|g (2) | < det(A) = E~ l M we will have 
Li ^> L 2 and L\L 2 ^> EM' 1 . We now write 

h® = E^(M^g? + s gf , gf) 

for % = 1,2. These vectors will then be a basis for Z 2 , and if x = Aih^+A 2 h^ 
is in the region given by \x\ — sox 2 \ < E X I 2 M~ X and \x 2 \ < E x l 2 then we 
will have |A,| < Li for i = 1,2. This allows us to make a change of basis, 
replacing (xi, x 2 ) by (Ai, A 2 ) so that our constraints on x\, x 2 are replaced by 
the conditions |A»| < L^. 

We may argue in exactly the same way for Z\ , z 2 using the fact that 

t = Zl /z 2 = 0(s Q ) + u4>'(s )u + O^x- 1 ) = 4>(s ) + 0(M~ l ). 

This allows us to replace the variables -21,-22 by Ti,t 2 subject to \t<\ < Tj. 
Here T\ ^> T 2 and T{F 2 <C FM~ X . These substitutions convert (j4j) into a 
new equation 

G (A 1 ,A 2 ;r 1 ,r 2 ) = l (11) 

say, where Go is bi-homogeneous of degree (2,1). Similarly they will turn 
(fTUj) into an equation of the shape 

G 1 (A 1 ,A 2 ;r 1 ,r 2 ) =0 (12) 

where G\ is bi-homogeneous of degree (a, b), say. Of course it is apparent 
from (TTTjl that the vectors (Ai, X 2 ) and (ti,t 2 ) will be primitive. 

When min(a, b) > 2 we can get a satisfactory bound from the following 
general result, which will be proved later, in £|6j 

Lemma 2 Let G{x\, x 2 ; yi, y 2 ) £ x 2 , yi, y 2 ] 6e an absolutely irreducible 
bi-homogeneous polynomial of degree (a, 6) a,b > 1. Le£ e > 6e given. 
Then for any X > 1 taere are O a .& )£ (X 2 / b+£ | |G| | £ ) points (xi, x 2 , yi, y 2 ) G Z 4 
satisfying the conditions 

g.c.d.(xi,x 2 ) = 1, g.c.d.(yi,y 2 ) = 1, 
G(xi, x 2 ; yi, y 2 ) = 0, and max|xj|<X. 
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Notice here that there is no size constraint on y± or y 2 . 

When a > 2 the lemma shows that ([T21 has 0^(T 1 1+,? a; T ') solutions Ai, A2, 
ri, r 2 . Each of these corresponds to at most one solution of (4), and therefore 
contributes 0(1) to M(E,F). Similarly, if b > 2 that there are O^L^x 11 ) 
solutions. 

We next dispose of the case in which a or b is zero. For example, if a = 
then (TT2"]) specifies a finite number O v (l) of pairs r 1; r 2 , and for each of these 
there is a corresponding pair (zi, z 2 ), producing a value of / = z\ + z\. Each 
such / contributes O(logo;) to Ai(E,F) by Estermann's bound ([3]). Thus 
the interval / contributes O v (logx) when a = 0. Similarly if 6 = there are 
0^(1) corresponding values for e. As in (jSJ) each such value e contributes 
<C p[e){xe~ 2 + 1) to Ai(E,F). This is also satisfactory, since 

e>E~>D = x 1/2 

and p(e) <C,, x v . Thus we have 0^(x' 7 ) solutions corresponding to I when 
min(a, b) = 0. 

When a = 1 the equation (1T21) can be written 

AiGn(ri, t 2 ) + X 2 G 12 {r 1 , t 2 ) = 0. 

Thus Ai = q^Guin, t 2 ), X 2 = -q~ 1 G n (T 1 ,T 2 ), where q divides Gn(n,T 2 ) 
and Gx 2 (ti,t 2 ). Since n and r 2 are coprime it follows that g divides the 
resolvent R of Gn and G\ 2 . This resolvent is non-zero since G\ is irreducible. 
Moreover it is bounded by a power of x, whence there are O v (x v ) possible 
choices for q. (The reader should recall at this point that the forms Gn and 
Gi 2 are determined, up to 0^(1) possibilities, by the interval J.) For each 
available choice of q we substitute our values for Ai, X 2 into ffTTj) to obtain a 
Thue equation G$(t\, t 2 ) = q 2 . Unfortunately we cannot use the full force of 
known results on such equations, since it is possible for G 3 to be a power of 
a linear form. None the less there can be at most 0(T{) possible pairs Tx,t 2 . 
It follows that we have at most O^x^Ti) solutions in total. The case 6 = 1 
is entirely analogous, leading to a bound O^x^Li). 

In summary we have a bound 0^(T 1 1+r? x' 7 ) on the number of solutions, in 
each of the cases a > 2, or a = 0, or a = 1. Similarly we have an estimate 
O v (L\ +t> x v ) whatever the value of b. We therefore conclude as follows. 

Lemma 3 For any 77 > the contribution to Ai(E,F) corresponding to a 
single interval I is O v (x 1] min(L} +,? , T 1 1+,? )). 
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5 Completion of the Proof 

Having fixed M as in Lemma 1 we must now sum up min(Li,Ti) for the 
various intervals /. In the notation of the previous section, if (0:1,22) corre- 
sponds to g (1) then L 1 (x 1 - s x 2 ) < M _1 a/E and L x x 2 < y/E. If L x > \/E 
we see that x 2 = 0, and then x\ = 0, which is impossible. The intervals 
/ = (s , s + M _1 ] will be produced by taking s = x 3 M~ l for integers 
£3 <C M. Thus the number of intervals for which L < L\ < 2L is at most 
the number of triples (xi,x 2 ,x 3 ) G Z 3 with g.c.d.(xi, x 2 ) = 1, for which 

x 2 x 3 = Mxi + 0(LT l y/E), x 2 < L^x/E, and x 3 < M. 

We now recall that L x > L 2 and that L X L 2 > EM' 1 . Thus L > £' 1 / 2 M _1 / 2 . 
Moreover, as noted above, we have L <C E 1 ! 2 . In particular, if M is large 
enough we can have x 2 x 3 = only when x\ = 0. Since X\ and x 2 are coprime 
this case can arise only when x 2 = ±1 and x 3 = 0. When x 2 x 3 7^ the 
conditions on x 2 and x 3 imply that X\ <C L^y/E, and a divisor function 
estimate then shows that there are O n (x ri L^ 1 y/E) pairs x 2 ,x 3 for each value 
of X\. We conclude that there are O r/ (x v L~ 2 E) intervals / for which L\ is of 
order L. Since each interval makes a contribution <C^ x v L\ +v } by Lemma [21 
we get a total <^ x^Lr^E < x 2n L' l E, since L < £ 1/2 < x. By dyadic 
subdivision for L > E 1 ' 2 ^!- 1 ' 2 we find that M(E,F) <^ x 2ri (EMfl 2 . 

We can prove a precisely analogous estimate JA(E,F) <C, x 2r] (FM) l l 2 
by considering the number of intervals J = (0(so), 0(so) + 0(M -1 )] which 
produce a value T x in a given dyadic range (T, 2T]. Here we use the fact 
that J C (t 3 M- 1 ,t 3 M~ 1 + 0(M -1 )] for some integer t 3 . We also need to 
remark that each value of £3 occurs 0(1) times, since |0'(s)| ^> 1 for the 
values of s under consideration. With these observations the argument then 
goes through just as before. We may therefore conclude that 



It remains to use this result with the value for M coming from Lemma 
CD It is convenient to write E = x^, so that F = x 2 ~ 2 ^ + °^ l / Xozx \ In view of 
our remarks at the end of §2 we have (essentially) 1/2 < ip < 3/4. We may 
then employ a value M with 



M{E, F) x 2v (mm(E, F)M) 



1/2 




This value will automatically satisfy M G [x 



/ 2 ,x] if 77 is small enough. It 



11 



follows that 

log M(E,F) ^ n 1 . .. n n/ , M . f 9-0(1 - -0) 1 

, < 2r]+ mm ^,2-2^ + (1 + 77 max — ^— — , - 

log a; ^ [84 

+ 0,((logx)- 1 ). 

However, since 

i m in W ,2-2^) + max{^»,i}<l 

for the relevant range of 0, we deduce that M.{E,F) <C, x 3??+7 / 12 , and our 
theorem then follows. 



6 Lemma 2 



Lemma [2] is closely related to two results of Broberg. In [2, Theorem 1] 
Broberg establishes a general result about finite covers of P 1 which, when 
translated into our notation, would provide an estimate OG, a ,b,e(X 2 / b+e ) of 
the desired order, but without any explicit dependence on G. This explicit 
dependence can be deduced from a second result of Broberg [3] , but this has 
not been formally published. We therefore give a brief sketch of a direct ar- 
gument independent of these two papers. This uses the determinant method, 
for more details of which the reader should consult [5], §3]. 
We will need a crude bound on the size of y. Let 



G(x;y) = G (x)y b 1 + ... + G b (x)y. 



The form Go cannot vanish identically since G is irreducible. Thus there are 
O a (l) primitive integer vectors x for which Go(x) = 0. It is not possible for 
all the forms Gj(x) to vanish simultaneously for a vector x ^ 0, since G is 
irreducible. Thus, with O aj b(l) exceptions, any solution of G(x;y) = has 
?/2 1 Go (x), with Go(x) ^ 0. We may therefore assume that 1 2/2 1 X a ||G||, 
and similarly for y\. It will be convenient to write these bounds in the form 
M, 1 2/2 1 < Y, with Y < X a ||G||. 

Our overall plan now is to apply the p-adic determinant method. By 
making an invertible integral linear substitution on x we may assume that 
Go(l, 0) 7^ 0. Indeed we can choose the coefficients of the substitution to be 
bounded in terms of a alone, so that we may still assume that |x| X. 
Suppose we have a parameter P ^> a ,b logX||G||. Then for any solution x, y 
of G(x; y) = with |x| <^ a X, |y| Y, we either have 

dG(x;j) 
X2V2 — k = 

oyi 
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or there is a prime p G (P, 2P] not dividing x 2 y2{dG / dy\) . The first case 
immediately give us an auxiliary bi-homogeneous form H (x; y) not divisible 
by G, at which our solution also vanishes. In the alternative case the point 
(xi/x2, J/i/j/2) lies above a smooth F p -point on the curve G(s, l;t, 1) = 0. 
We can then expand 2/1/2/2 as a p-adic power series in xi/x 2 , as in Lemma 5 
of the author's paper [5]. We then consider the matrix of bi-homogeneous 
monomials in x and y of degree (H, b — 1). There are k := (H + 1)6 such 
monomials. The corresponding k x k determinant then has archimedean size 
^a,b,H (X H Y b ~ 1 ) k . Moreover it will be divisible by p k ( k ^ 1 )/ 2 . The argument 
of 0, §3] then produces an auxiliary form if(x;y) providing that 

We now recall that Y <C X a ||G||. Thus on choosing H sufficiently large 
we see that it suffices to have p ~^> a ,b,e (X\ \G\ \) £ X 2 l b . In addition to the 
form x 2 y 2 dG j dy\ that we have already mentioned we now obtain one further 
form H(x; y) for each F p -point on the curve G(s, 1; t, 1) = 0. We thus con- 
clude that every solution to G(x; y) = with max|xj| < X satisfies one of 
O a ^ £ {{X\\G\\) £ X 2 / b ) auxiliary conditions F(x;y) = 0. Here H is a bilinear 
form coprime to G, with degrees bounded in terms of a, b and e. For each such 
form, there are O a ,b,eiX) common solutions to G(s, 1; t, 1) = H(s, 1; £, 1) = 
by Bezout's Theorem, and s, t determine x, y since these vectors are primitive. 
This suffices for the lemma. 

7 Further Improvements 

It is possible to reduce slightly the exponent 7/12 occurring in the theorem. 
Since the improvement is very small we content ourselves with a very brief 
sketch of the argument. 

The first step is to repeat the analysis of §3 taking K = L = 1 and 
obtaining a bi-linear form F in (TTOj) . providing that M £ [x 1 ^ 2 ^] satisfies 
M > (£i/3_pi/2)i+<5_ Note here that in fact 

2.1/2 < (£1/3^1/2)1+5 < x 

for large enough x and small enough S, since x 2 <§C E 2 F x 2 and E,F ^> 

When F is bi-linear the Thue equation referred to in §4 will have degree 
3, and will produce O v (x v ) solutions for each interval /, except when G3 is 
proportional to a cube. In this case the corresponding solutions of 

2xxx 2 yi + {x\ - x 2 2 )y 2 = 1 (13) 
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lie on a line 

( Xl ,x 2 , yi ,y 2 ) = {xf\xf\yf\y^) + X^,^ M ) 



contained in the variety f|T3|) . 

We now write A4o(E, F) for the number of quadruples {x\,x%, yi, y 2 ) sat- 
isfying f|T3|) . for which 

E/2<xl + x 2 2 <E, F/2<yl + y 2 2 <F 

but which do not lie on a line in the variety f|T3|) . We may then deduce that 



M {E,F) <„ x^E^F 1 ' 2 < <r? x 2 " +1 - 2 ^ 3 , 

where E = as before. Alternatively we can use our previous argument 
which shows that 

M {E,F) < M(E, F) <„ x 3^+min(^,2-2^)/2+9^(l~^)/8_ 

These suffice to show that 

M (E,F) <£ v x 3r,+ ™ 

with 

26W433 = 577g8 (< 1 = 0.58333...) 
81 v 12 y 



the critical value of ^ being (55 — V433)/54. 

It then remains to consider the form taken by lines lying in the surface 
(|T3|) . The lines which contain more than one integral point may be described 
explicitly, and one is then able to show that they contribute O v (x ri+1 ^ 2 ) to 
M(E,F). 

In this way one may improve the exponent in the theorem to w. 
The author is grateful to Thomas Reuss for pointing out an error in the 
original version of this final section. 
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